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ABSTRACT 

Analysis of galaxies with overlapping images offers a direct way to probe the distribution 
of dust extinction and its effects on the background light. We present a catalog of 1990 such 
galaxy pairs selected from the Sloan Digital Sky Survey (SDSS) by volunteers of the Galaxy Zoo 
project. We highlight subsamples which are particularly useful for retrieving such properties of 
the dust distribution as UV extinction, the extent perpendicular to the disk plane, and extinction 
in the inner parts of disks. The sample spans wide ranges of morphology and surface brightness, 
opening up the possibility of using this technique to address systematic changes in dust extinction 
or distribution with galaxy type. This sample will form the basis for forthcoming work on the 
ranges of dust distributions in local disk galaxies, both for their astrophysical implications and as 
the low-redshift part of a study of the evolution of dust properties. Separate lists and figures show 
deep overlaps, where the inner regions of the foreground galaxy are backlit, and the relatively 
small number of previously- known overlapping pairs outside the SDSS DR7 sky coverage. 



Subject headings: galaxies: spiral — dust — galaxies: ISM 



1. Introduction 

Dust composes only a small fraction of the 
mass of the interstellar medium in galaxies, but 
it plays key roles in both their development and 
observed properties. Radiative cooling by dust 
facilitates the collapse of massive clouds during 
episodes of star formation, and grains remain as 
a key constituent of protoplanetary disks around 
the nascent stars. At the other end of the stars 
lives, grains condense in the cooling and enriched 
gaseous envelopes of supergiants, planetary nebu- 
lae, and the eject a of novae and supernovae. On 
grander scales, dust influences the evolution of 
galaxies even as its absorption complicates our 
ability to interpret the observed properties of 
galaxies. Correction for dust extinction within 
galaxies is important for the use of SN la as cos- 
mological probes, especially if the extinction curve 
varies among galaxies. It also affects distance 
measures via the Tully-Fisher relation, through 
the required corrections for internal extinction as 
a function of inclination to our line of sight. If 
the extinction is high enough toward the cores of 
luminous spirals, many straightforward estimates 
of stellar masses and distributions come into ques- 
tion ( Driver et al.||2QQ8 ). We present an extensive 
catalog of superimposed galaxy pairs, which are 
useful for measuring dust extinction and allow a 
much wider study than previously available. 

For decades, the standard view on dust effects 



scale height of the dust is much less than that 
of the stars. Around the same time, |Valentijn| 



in galaxies derived from variants of the Holmberg 



( 1958 ) approach, based on surface brightnesses of 
disks viewed at different inclinations to the line of 
sight. Optically thin disks will increase in surface 
brightness when viewed closer to edge-on, while 
in the optically thick limit, the surface brightness 
will remain constant. Holmberg's data indicated a 
mild inclination dependence and correspondingly 
modest internal extinction. This comforting con- 
clusion was eventually challenged by [Disney et al 
(1989) on theoretical grounds; models with very 



dusty disks could reproduce the data if the disk 
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(1990) claimed that extensive analysis of surface 
photometry of galaxies from the ESO-LV sample 
suggested nearly constant surface brightness with 
disk inclination, and much higher internal opacity 
than previously assumed. These results spurred a 
rebirth in work on this question, leading to a wide 
range of studies converging on a general picture 
of high opacity in spiral arms, resonance rings, 
and the centers of disks, with diffuse dust fading 
to very small extinction at the edges of the opti- 
cal disks (Thronson 1995). Still, it became clear 



that substantial variations exist among and within 
galaxies. Use of the Galaxy Zoo morphological 
classifications with color data suggested further 
complications, such as very luminous disks being 
comparatively dust- deficient ( Masters et al.|2Ql'Q ). 

All these factors can be exacerbated for galax- 
ies at high redshift. Many have high rates of star 
formation, suggesting a rich and massive interstel- 
lar medium, and are observed most often in the 
emitted ultraviolet, where dust effects are more 
important, compared to the emitted optical range 
where our basis for comparison is more extensive. 
Different ways of correcting for internal extinction 
give widely varying results, contributing to the dis- 
parate conclusions obtained by several groups on 
the history of cosmic star formation, even when 
starting from the same Hubble data in the deep 



fields (e.g., Madau et al. 


1996, 


Thompson, Wey- 


mann, & Storrie-Lombardi||2001 


). 



The dust content and distribution in galaxies 
may be approached in several ways. Together with 
models of radiative transfer and constraints from 
the global energy budget, fits to the spectral en- 
ergy distributions (SEDs) of galaxies have been 
used to estimate the typical extinction at various 
wavelengths and total dust masses. Models must 
assume the relative distributions of stars and dust, 
as set out by, for example, Witt, Thronson, fc 



Capuano (1992); dust mixed with the stars gives 
much less overall extinction than a simple fore- 
ground screen. This must be taken into account 
in deriving extinction laws from comparison of ob- 
jects' spectra, as in the effective extinction law 
derived by'Calzetti, Kinney, fc Storchi-Bergmann| 
( |l999) . SED fitting of edge-on spirals consistently 
underestimates the face-on extinction and amount 



of UV radiation which is reprocessed (Baes et al. 
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2003 Baes et al. 2010, Kuchinski et al. 1998), so 



improved measurements of these quantities are im- 
portant in refining models of this kind. 

Far-infrared emission provides very direct map- 
ping of the dust content of galaxies, modulated by 
the typical grain temperature. Very cold grains 
remain difficult to detect, and large fractions of 
the total dust mass can in principle be found only 
by careful spectral decomposition of the submil- 
limetre emission, in view of the dominant contri- 
bution of the warmest grains. Far-IR data have 
hinted that some spirals have very extensive dust 
distributions, at least through a larger disk scale 
length for dust than for starlight (Davies et al. 
[1999, S un fc Hirashita 2011) . Aside from the 
detectability of dust being strongly weighted to 
warmer grain populations, far-IR emission studies 
still suffer from having spatial resolutions much 
poorer than is available using optical techniques. 

Herschel observations have contributed consid- 
erably to understanding the radial and temper- 
ature distributions of galactic dust, and the to- 
tal mass in grains. With adequate spatial reso- 
lution, the mix of grain temperatures should be 
coupled to location in the galaxy through the ra- 
diation field of starlight. Herschel traces dust at 
350 /im out to 0.8-1.3 times the diameter D25 of 
the de Vaucouleurs isophote in luminous spirals 
dBendo et al.|2010[ [Pohlen et al.pOlOl ). While the 
gas/dust ratio increases with radius in Virgo spi- 
rals, as might be expected from metallicity gra- 
dients (jSmith et all 120101 iPohlen et al.l 120101), 



across whole galaxies (Alton et al. 2000, James 
eTaLll 2002D . 



the dust shows a cutoff matching that of H I in 
stripped galaxies (Cortese et al. 2010), testifying 



to a dynamical link between gas and dust. Typical 
Herschel observations reach dust column densities 
close to 0.05 Mq pc''^ at T = 19 - 22 K (e.g. 
Smith et aLpOK T). Even with these dramatic im- 



provements in sensitivity and angular resolution, 
degeneracies remain in interpreting the grain pop- 
ulations in galaxies. These arise largely because 
changes in the emissivity parameter P can mimic 
changes in the relative contributions of warm and 
cold dust across the FIR and submm regimes. The 
value of P expresses departures from the black- 
body intensity Bx{T) for an emitter which is small 
compared to the relevant wavelengths, of the form 
/ oc X~^Bx{T). Changes in the distribution of 
grain size will affect /3, although it remains unclear 
whether this would be measurable when averaged 



Absorption studies have the high angular res- 
olution of optical or ultraviolet imaging, but are 
limited to regions with background sources that 
are adequate in solid- angle coverage and our un- 
derstanding of their properties. When the back- 
ground light arises within the galaxy under study, 
other issues arise in our knowledge of the rela- 
tive distribution of stars and dust, and the role of 
scattering; these are approached in such ways as 
modeling the ^-distribution of disk stars of various 
kinds to interpret multicolor absorption measure- 
ments (Elmegreen 1980). In comparison, use of 



a more distant galaxy as the backlighting source 
gives several advantages, at the expense of dra- 
matically reducing the range of galaxies available 
for analysis. The spatial resolution is limited only 
by the telescope's image quality, and we need know 
nothing about the internal stellar structure of the 
foreground galaxy. Scattering corrections are neg- 
ligible once the two galaxies have a line-of-sight 
separation only a few times their diameters, sat- 
isfied if their redshift difference indicates that the 
galaxies are far apart along the line of sight. High- 
quality extinction measurements can be obtained 
even in the outskirts of a galaxy disk. 

Results of both these approaches to absorption 
will have a clumping dependence; internal struc- 
ture in absorption can make small clumps harder 
to detect if they are below the spatial resolution of 
the data, and the derived reddening law will gen- 
erally be ffatter (grayer) than the intrinsic form 
given by the grain properties, due to mixing of 
regions with different transparency within a sin- 
gle resolution element. The dust clumping also 
enters into the distinction between dust mean col- 
umn density as derived from modeling FIR mea- 
surements and the optical extinction we measure 
here. 

Application of the overlap technique has im- 
proved with data quality. |Keel| ( |1983| ) showed a 
first application to NGC 3314 and the foreground 
system in NGC 1275, with vidicon imagery of lim- 
ited signal-to-noise ratio. This approach was put 
on a genuinely quantitative footing with analy- 
ses of CCD imagery by [Andredakis fc van der 
( [T992I), [White fc Keel| ( |1992D,|Berlind et al. 
(|1997|, [White et al.| (|2000D and [Domingue et al. 



(1999) (D1999) , in the last case comparing optical 
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extinction with far-IR and submillimeter data to 
compare dust masses estimated from emission and 
absorption. Broadly, these studies show that ex- 
tinction may be high {Ab > 1) within spiral arms 
at a wide range of galactocentric radii, with inter- 
arm extinction smaller and declining outward in 
a roughly exponential manner. Extinction in the 
inner kpc has been measured only in NGC 3314a, 
where it reaches > 5. One resonance-ring spi- 
ral was observed, in which the ring opacity is much 
larger than that found on either side of it. Extend- 
ing the analysis to the resolution of tens of par- 



sees enabled by the Hubble Space Telescope (Keel 



&; WhitepQOTal |Keel White||2QQlb| [Elmegreen 



et al.||2QQl[ [Holwerda et al.| |2QQ9D showed that 
the effective reddening law depends significantly 
on linear resolution. For a clumpy distribution 
of dust, this is unavoidable, since the weighting 
of regions in transmission becomes wavelength- 
dependent ( Fischera, Dopita, fc Sutherland|2Q03 ). 
The values measured in the outer disks of sev- 
eral spirals at HST resolution approach the local 
Milky Way mean, suggesting that the extinction 
contrast of dust structure on these scales is mod- 
est enough that we may plausibly be approaching 
the true grain reddening behavior. These studies 
were all limited by the small number of suitable 
and nearby galaxy pairs known, so that the be- 
havior with galaxy type and its variance could not 
be explored. In particular, only the most sym- 
metric spirals - grand-design and strongly ringed 
systems - could be measured when the sample size 
was too small to average over unmodelled struc- 
ture in individual galaxy disks. 

The ideal pair for mapping dust would consist 
of a face-on, fairly symmetric spiral, seen nearly 
half backlit by a smooth elliptical or SO system, 
of which half is seen free of any foreground ex- 
tinction. In such a system, there is a large region 
in which each galaxy is seen essentially by itself 
so that a good model of each galaxy can be pro- 
duced, and a large region over which dust can be 
mapped. If B is the background intensity and F 
the foreground intensity, modelled point-by-point, 
the optical depth r may be estimated from the 
observed intensity / using 

--^ 

Various compromises may have to be allowed for 
some uses; for example, only in spiral/spiral over- 



laps is there enough fiux to measure extinction 
into the ultraviolet, although the point-by-point 
errors are unavoidably large due to the rich struc- 
ture of spirals which is most pronounced at such 
short wavelengths. If they are seen through par- 
ticular parts of a spiral, even background galaxies 
of comparatively small angular size may be useful 
(if seen at the disk edge or across a spiral arm, for 
example); selection effects dictate care in analy- 
sis because of the small objects that will not be 



seen through relatively opaque regions (Holwerda 



et al. 2007c). Useful spectroscopic estimates of 



extinction can be made even when there are sig- 
nificant departures from symmetry, using the rel- 
ative amounts of light at each galaxy's redshift 
dPomingue et al.||20Q0[ D2000). 

The ability to harvest large samples of over- 
lapping galaxies from surveys like the SDSS has 
reinvigorated the study of extinction using back- 
lighting. We present here a new, large listing of 
galaxy pairs suitable for such dust studies. Its 
production has relied on contributions from many 
of the volunteer participants in the Galaxy Zoo 
project ( Lintott et al. 2008) . We will use this sam- 
ple to address a range of issues in galaxy dust con- 
tent via this single technique. What are the sys- 
tematics with Hubble type, central surface bright- 
ness, strength of bars or rings? How common are 
the kinds of extended dusty disks found by |Hol- 



werda et al. (2009) ? Is there a very dusty mor- 
phological disk sequence which doesn't stand out 
in color alone ("Thronson 1995)? In addition to 



the study of dust extinction in a variety of envi- 
ronments, uses of the listing range from testing 
image-decomposition routines to correlation with 
supernova detections as a way of measuring red- 
dening laws. 

2. Galaxy Selection 

The bulk of this catalog consists of pairs origi- 
nally selected by volunteer member s of the Galaxy 
Zoo project^ As described by Lintott et al. 
( [2Q08| ), the main goals of the project are reached 
by visual classification of galaxies from the SDSS 
via a web interface. This project operated outside 
the primary statistically-oriented framework of the 
Galaxy Zoo project, as so-called forum science. 
After noticing that some of the pairs being shown 



■"^http:/ /www. galaxyzoo.org 
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and discussed on the project forunj^ were likely 
non-interacting, overlapping pairs, WCK posted 
a specific request for such systems on the forum. 
Users responded enthusiastically, posting images 
and identifications for large numbers of candi- 
date^ We inspected these candidate pairs (and 
objects appearing on other discussion threads of 
the forum, a total of ~ 7000 suggested pairs) for 
overlap, utility for dust study, and evidence of in- 
teraction or other asymmetry, with surviving can- 
didates going into our working list. Galaxy Zoo vi- 
sual detections are now by far the dominant source 
for low-redshift overlapping systems. Prior to sur- 
veys from the SDSS images dHolwerda et al.|2007b 
and this work), no more than 25 such pairs were 
mentioned in the literature (Appendix A). 

For completeness, we checked for pairs already 
selected spectroscopically from the SDSS (Holw- 



erda et al. 2007b). Modifying the criteria used 



for the gravitational-lens search by Bolt on et al 
(2004), that project selected objects where the 
SDSS spectra showed an emission-line redshift 
lower than the cross-correlation (absorption-line) 
value. This target list was winnowed visually 
for appropriate geometry (images actually over- 
lapping) to give a final set of 83 pairs. Of these, 
34 are not useful morphologically for our sample 
here, and 9 were already selected by the Galaxy 
Zoo project. The remainder were objects in which 
two galaxies are not apparent on inspection of 
the SDSS images, such as having a low-luminosity 
foreground object closely aligned with a luminous 
background galaxy (so our image analysis is not 
appropriate, at least with ground-based data). 

Some of the newly-recognized systems are large 
enough and bright enough to appear in earlier sur- 
veys (such as the Palomar Sky Survey), but the 
depth of the SDSS data and the helpful genera- 
tion of color-composite images with high-dynamic- 
range mapping (| Lupton et al.|[2QQ4 ) made it easy 
to recognize features that were overexposed on 
most digital representations of the Palomar pho- 
tographs. The color information was especially 
valuable for finding reddened galaxies behind the 
disks of spirals, which may be prominent only in 
the i and z bands. 



^http:/ /www. galaxyzooforum.org 

^The dedicated discussion thread for these pairs is 
jhttp: / / www.galaxyzooforum.org/index.php?topic=6732.0| 



As an additional check for automated selec- 
tion of overlapping pairs. Galaxy Zoo volunteer 
Lionel Poffel contributed the results of an SQL 
query to identify pairs of objects closer than the 
sum of their listed Petrosian radii, with redshift 
difference Az > 0.01. With hindsight, this ap- 
proach was rendered less effective by blending is- 
sues; for strongly overlapping galaxy images, the 
SDSS photometric pipeline often extracts only the 
nucleus of one galaxy, assigning a Petrosian radius 
too small by factors 4-5. A converse issue is that 
stars can be misclassified as galaxies when super- 
imposed on the outskirts of a large galaxy image. 
Indeed, of 120 candidate pairs from this selection, 
24 involve a foreground star. Even after remov- 
ing these false alarms, this query added 28 new 
pairs to the catalog, and recovered 14 otherwise 
selected (plus 49 objects which had been consid- 
ered from forum postings and rejected for reasons 
of symmetry, low surface brightness in the overlap- 
ping regions, or inappropriate orientations). This 
comparison testifies to the utility of visual selec- 
tion by the Galaxy Zoo volunteers in generating a 
candidate sample for this project. 

Not all galaxies in our catalog are in the SDSS 
database. A few were otherwise known from ear- 
lier results and lie outside its footprint (Appendix 
A). Even within the imaging coverage, some have 
at least one member not identified as a separate 
photometric object. 

Of 1990 pairs comprising our final catalog, 1713 
(86%) have one or more known redshifts mea- 
sured in the SDSS or tabulated by NEE|5 The 
NED additions are for bright galaxies or com- 
panions to galaxies with SDSS spectra, where 
there is no SDSS spectrum because of sampling 
or fiber-collision limitations. We adopt a redshift- 
difference criterion to fiat the subset of pairs which 
cannot be physically interacting, driven by the 
range of Az seen in tidally-interacting galaxy pairs 
and including high- velocity encounters in clusters. 
We adopt Az > 0.008 {Av > 2400 km s"\ a very 
conservative value which would distinguish non- 
associated galaxies even in the presence of signif- 
icant redshift errors. Of 494 pairs with both red- 
shifts known, 218 satisfy this criterion, ruling out 
gravitational interaction with each other. For this 
subset, scattering effects are guaranteed to be neg- 



^ http://ned.ipac. caltech.edu/ 
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ligible in measurements of obscuration ( White et 
al.|2QQQ ), and there wih be no asymmetries driven 
by interaction between the pair members. We re- 
jected candidate pairs which have these large red- 
shift differences, but in the wrong sense to be use- 
ful for dust studies (such as an elliptical galaxy in 
the foreground of a spiral). 

The catalog properties in distance and Az are 
illustrated in Fig. [l] which shows a histogram 
of the smaller redshift value Zmin^ when both are 
known, or the single known redshift, and plots val- 
ues when both Zmin and Zmax are known. For sys- 
tems with one redshift known, the median value 
is z = 0.070, with quartiles 0.041, 0.105. Combin- 
ing single redshifts with foreground values when 
both are known, the median is z = 0.064 (quar- 
tiles 0.036, 0.097), while for background systems 
when both values are known, the range rises as 
expected: median z = 0.077 and quartiles 0.040, 
0.127. 

WCK screened pairs posted by Galaxy Zoo 
participants for suitability, in particular rejecting 
pairs with obvious tidal interaction as seen on the 
SDSS color composite images. In a final round, 
WCK and AMM independently re-examined can- 
didates, to check for clearly interacting pairs which 
had survived the initial round. At this stage, we 
rejected about 10% of systems initially selected for 
the catalog. These were pairs with obvious tidal 
tails or decentered nuclei, objects where it was not 
clear whether there are in fact two distinct galaxies 
rather than unusual substructure in a single one, 
and pairs in which the redshift difference is both in 
the wrong sense to be useful for dust studies (spi- 
ral behind elliptical) and large enough to clearly 
show that the galaxies are not in the same group, 
so that the redshifts will indicate distance order- 
ing. This "discard" list is available in case some 
of these systems prove of use for related questions. 
We required positive evidence of tidal distortion to 
reject a candidate pair; some interacting systems 
with only weak distortion, or tidal features evident 
only at low surface brightness, will remain in the 
sample. 

Our catalog of 1990 overlapping galaxy pairs 
is listed in Table [l] The entire listing is given 
in the online edition; the printed version includes 
a subset to illustrate content and format. This 
subset includes the initial lines of the RA-ordered 
list, all the newly- found objects shown in Fig. 
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Fig. 1. — Distribution of known redshifts of cat- 
alogued pair members. Points indicate lower and 
higher redshift values Zmin and z^ax when both 
are known. In the upper panel, the superimposed 
histogram shows the distribution of single redshift 
values and Zmin in bins of Az=0.02, with values 
scaled down by factor 500 to fit the coordinate 
axes. The tight grouping of points along the di- 
agonal marks pairs with close enough redshifts to 
be physically associated on group scales, although 
not necessarily interacting. This is emphasized in 
the lower panel, which shows the individual pair 
differences in radial velocity Av. The dotted hori- 
zontal line marks our Av = 2400 km s~^ division; 
pairs above this are not physically associated. 
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[2j and a few additional objects to illustrate the 
combinations of data and identifications found. 
Data for previously known pairs in Fig. [2] are 
listed in Table [4j SDSS pairs are listed with 
the coordinate designation of the brighter mem- 
ber. Magnitudes are SDSS inodelinag_r values in 
the r band. These must be treated with cau- 
tion, if not suspicion, because pipeline separation 
of strongly blended galaxy images often results in 
assigning much of the fainter galaxy's flux to the 
brighter one (section 2.2), and spot checks also 
show some pairs in which even a smooth symmet- 
ric galaxy is broken into multiple SDSS photomet- 
ric objects.. A second magnitude with only one 
significant figure past the decimal was estimated 
from the SDSS composite image, if the galaxy was 
not detected as an SDSS object (or as multiple 
objects). A value of 30.00 was assigned if the im- 
ages are so strongly blended that no estimate for 
the fainter galaxy was reasonable. In some cases 
there are multiple background galaxies; the mag- 
nitude listed is for the brightest of these. We in- 
clude redshifts where features for both members 
were clearly detected in a single SDSS spectrum 
(indicated by fg and bg for fore/background). The 
catalog in sortable web form, with additional PDF 
files with finding charts and additional data on 
each pair as well as our "reject" list, is available 



from http: / / data.zooniverse.org/ overlaps.html, 



2.1. Pair categories 

We note some special categories defined by ge- 
ometry or redshifts with a "type" indicator in Ta- 
ble [l] which mark subsets of the catalog particu- 
larly useful for various purposes. They are denoted 
by mnemonic designations, which are summarized 
along with their frequency in the catalog in Table 
[2] A selection of each is illustrated in Fig. [2] 

F: spirals seen nearly face-on in front of an el- 
liptical or SO background system. These are clos- 
est to the ideal for most kinds of extinction stud- 
ies. For our purposes, these include spirals face-on 
enough that the extinction structure is dominated 
by arm/interarm variations, rather than the ex- 
tent of dust perpendicular to the disk. In practice 
this includes spirals with planes inclined as much 
as 60° to the plane of the sky. 

Q: the background galaxy is nearly edge-on 
and is projected nearly radial to the foreground 
galaxy, so the backlit area spans a large range in 



radius but a narrow one azimuthally within the 
foreground system. This means that symmetry 
requirements on the foreground galaxy can be re- 
laxed substantially, modeling it by interpolation 
across a narrow angular sector. Arp 198 (UGC 
6073, VV 267; Fig. |2| is a good example, in 
which the run of extinction with radius can be 
retrieved even in the presence of a rather com- 
plex foreground spiral pattern. In the most fa- 
vorable cases, the SDSS images show evidence of 
extinction to indicate that the edge-on galaxy is in 
the background; better imagery may show some of 
the fainter Q systems to in fact have the edge-on 
galaxy in front. 

^: the spiral is seen essentially edge-on, at least 
partially backlit by a smooth galaxy. These are 
useful for studying the structure of extinction per- 
pendicular to the disk of the foreground galaxy, 
and the radial extent of extinction to very sensi- 
tive levels. The letter is selected to remind one of 
a thin disk with a round galaxy behind it. 

X: both galaxies are seen nearly edge-on, with 
their disks crossing as seen either near both nuclei 
or along one disk. These are relevant to the occur- 
rence of extended ^-distributions of dust (as in ^ 
types), as well as to questions of the distributions 
of angular momenta in pairs (with geometrically- 
based corrections for the likelihood of seeing the 
disks intersect when they have various actual an- 
gular differences). A handful of these (coded with 
R, as below) may in fact be polar-ring systems 
seen edge-on to both the central galaxy and the 
ring. 

SE: spiral/elliptical superpositions that do not 
fall in one of the other geometric categories above. 

S: spiral/spiral overlaps. These have much 
richer background structure than S/E overlaps, 
but are useful in probing extinction into the ultra- 
violet because the background source will remain 
detectable at much shorter wavelengths. At least 
using averaging techniques, they can be used to 
compare optical and UV extinction. With high- 
quality images, additional symmetry clues such 
as the direction of arm features can help sepa- 



rate foreground and background structure (Keel 



& White 2001b). 



B: the background galaxy has much smaller an- 
gular size than the foreground disk. Such align- 
ments are universal for small enough background 
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systems (transitioning, when their statistics are 
large, into the "synthetic field" technique match- 
ing galaxy counts and colors to global means and 
taking into account the cosmic variance; 



Holwerda 



etlir][2QQ5al [Holwerda et aL||2QQ5b| [Holwerda et 
al.|2QQ5c| " Holwerda et al.|2QQ7a ). As these papers 



demonstrate, considerable care is needed using nu- 
merous small background galaxies to map extinc- 
tion, because of bias in favour of the more trans- 
parent regions. Therefore, our listing includes only 
those where the line of sight to the background 
galaxy is particularly useful as a dust probe - edge 
of the optical disk, inside a resonance ring, or be- 
hind an outer spiral arm. In these cases even a 
local measurement of differential extinction can 
be useful. This category in particular employed 
subjective criteria as to utility. 

E: pairs containing only elliptical or SO galaxies, 
as judged from the SDSS. Our search concentrated 
on spirals, so coverage of these systems is much 
less complete. Many such pairs occur in clusters, 
but the possibility of subtle tidal distortions makes 
use of pairs with large redshift differences more 
valuable in seeking even low-opacity diffuse dust 
components. 

R: we use this code to flag a few objects which 
might be either near-central superpositions or po- 
lar rings, where better imagery or spectroscopic 
information would be needed to be sure. All of 
these, if overlaps, have very small central impact 
parameters and also appear in our list of deep over- 
laps. 

Nineteen pairs do not fall in any of these 
geometrically-defined categories, and are labelled 
simply "misc" in the table. 

Az: these are selected to have a known red- 
shift difference so large that the two galaxies will 
not be interacting with each other. Starting with 
known pairs, we take Az = 0.008 as the demarca- 
tion point. These pairs are the least likely to show 
tidal features or have tidally-induced asymmetries, 
and are so far apart that scattering corrections are 
negligible ( White et al. 2000 ). This category is 
independent of the morphological categories, and 
subsets of each type fulfill the Az criterion. The 
electronic table lists this in a separate column for 
ease of search. 

Among several of these categories we find deep 
overlaps, pairs in which the background galaxy is 



seen through the inner part of the foreground sys- 
tem. The reddening is often so strong that they 
can be selected from the SDSS gri composites for 
showing the background nucleus only at i. The 



nearby type example is NGC 3314 (Keel & White 
[2001b ); this catalog includes several additional ex- 
amples of central spiral/spiral overlaps, as well as 
some fairly nearby instances of backlighting ex- 
tending almost to the foreground nucleus. These 
objects are important (although we do not know 
how complete the selection is) because of evidence 
that the inner parts of some disks have significant 
optical depth even in the K band, making further 
IR observations a promising approach for extinc- 
tion measurements very deep inside these spirals. 
Table [3] lists such deep overlaps from our cata- 
log which are like NGC 3314 in the sense of hav- 
ing a background galaxy large and bright enough, 
projected close enough, for high-surface-brightness 
backlighting of the inner regions (typically inner 
3 kpc) of a foreground spiral which is not nearly 
edge-on. Six of these fall in our type R, which 
may also be polar ring systems. Table [3] includes 
the projected angular separation between galaxy 
nuclei. 

It is difficult to make an external assessment 
of the completeness of our catalog in magnitude, 
since the SDSS photometric pipeline sometimes 
fails in apportioning flux between deeply blended 
objects. Likewise, the pipeline size measurements, 
such as Petrosian radius, are not reliable for much 
of the sample. 

2.2. Redshifts 

Since photometric redshifts can be accurate 
enough, in some applications, to pick out galaxy 



pairs with matching redshifts ( Quadri & Williams 
2010), we have explored the use of the SDSS pho- 



tometric redshift estimators to distinguish addi- 
tional pairs with large redshift differences, with 
disappointing (although not wholly unexpected) 
results. For 156 distinct pairs having two spec- 
troscopic redshifts from SDSS DR7, we consid- 
ered the template-based and neural-network esti- 
mates photoz, photoz2 and their error estimates 
from the SDSS database. Of these, 33 showed 
catastrophic failures in both estimators, with er- 
ror > 5(7, traceable to deblending problems. In 
these cases, at least one pair member has its Pet- 
rosian radius underestimated by factors > 4, due 
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Fig. 2. — Sample pairs of various geometrical types. Images are from the JPEG files delivered by the SDSS 
SkyServer or our WIYN imaging displayed with simiS^r intensity mapping. North is at the top; image sizes 
are 1, 2, or 4 arcminutes square. Pairs are identified by common name or truncated coordinate designation 
to facilitate finding them in Table fl] (Continued on next page) 
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Fig. 3. — Deep overlaps, similar to NGC 3314 in having the galaxies' nuclei projected well within the area 
of the foreground disk, so there is a bright background source to study the extinction deep within a disk 
which is not itself edge-on (listed in Table [3|. Images are gri color composites from the SDSS Sky Server; 
all are 50" square except NGC 5021 and 6175 which are 100" square. Six of these fall in catalog category R, 
and might be similar to polar-ring systems rather than being two independent galaxies; redshift information 
would distinguish these. NGC 3314 itself, lying outside the SDSS region, is not shown here; it appears in 
Fig. Ill 
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to the deblending algorithm assigning only the nu- 
cleus to this object and the additional diffuse flux 
to the other galaxy or to neither. This generally 
carries corresponding photometric errors; a strong 
correlation between errors for the two estimators 
indicates that the photometric errors are mostly 
responsible. Even for pair members without ob- 
vious photometric problems, the photometric es- 
timates provide almost no information for z < 0.1 
(where both are biased upward, by 0.03 for pho- 
toz and 0.04 for photoz2)^ and are useful only for 
Az > 0.15 (in which case the apparent sizes of the 
galaxies otherwise suggest large differences inde- 
pendently). 

Our best assessment of the completeness of 
sample selection at this point comes from the 
subset with spectroscopic redshifts, which we can 
compare to the entire SDSS to ask what fraction 
of galaxies are selected as a function of z. Fig. [4] 
shows the fraction of SDSS galaxies in our cata- 
log as a function of redshift. We take the lower 
redshift of a pair when both values are known, 
since the foreground galaxy is the one under study. 
A galaxy is 10 times more likely to be selected 
at z = 0.01 as at z = 0.1, although the proba- 
bility stay roughly constant from z = 0.1 — 0.2. 
We can investigate how much of this effect is due 
to B-type pairs, with a faint background system; 
for more distant foreground systems, there is less 
volume available for background systems within 
the effective distance limit for seeing structure in 
SDSS images. This effect is significant only for 
z < 0.02; the lower trace in Fig. [4]shows the result 
of omitting the B-type foreground systems where 
the smaller system would be missed at greater dis- 
tances. This accounts for only a small fraction of 
the selection change with z. A large part of the 
effect must be in resolution, as the characteris- 
tic structures picking out two galaxies, and distin- 
guishing spirals from ellipticals, become so blurred 
that only more obvious pairings are found. 

2.3. Magnitude selection 

Statistics with pair and component magnitude 
provide some insight into the sample selection, 
subject to the caveats above about pipeline mag- 



is more robust to how the light was deblended 
into constituent galaxies. The magnitude differ- 
ence between pair members shows, as expected, 
a systematic decrease for fainter pairs, albeit with 
large scatter, the whole range r2— ri < 3.5 remains 
well populated for objects brighter than rtot = 18 
(Fig. H 

The counts of sample pairs with magnitude fol- 
low the Euclidean slope only at the bright end 
(Fig. [6|. When separated into bins by Am = 
^2 — ^1, the distributions remain broadly similar; 
incompleteness becomes serious at levels rtot = 
13 — 14, reaching fainter for more nearly equal 
pairs. Again, as one might expect, more equal 
pairs are selected preferentially among fainter sys- 
tems. 

3. Summary 

From thousands of candidates identified by 
Galaxy Zoo participants using the Sloan Digital 
Sky Survey, we have produced a catalog of 1990 
galaxy pairs with geometry of overlap suitable for 
use in the study of dust extinction. These are 
coded to distinguish various kinds of overlap use- 
ful for different aspects of the problem - spirals 
partially backlit by E or SO galaxies, spiral/spiral 
pairs, systems with the background galaxy edge- 
on and radial to the foreground system, and so 
on. About 11% of these have known redshift dif- 
ferences large enough to eliminate the possibil- 
ity of interaction with one another, and thus the 
possibility of such an interaction causing a break- 
down of the galaxy symmetry needed to analyze 
the extinction. Use of these pairs to study ex- 
tinction complements infrared techniques, giving 
the high angular resolution of optical observations 
and retaining sensitivity to even very cold grain 
populations. This sample can also be used to 
search for additional examples of very extended 
dust structures, such as found by |Holwerda et al.| 
(2009). Additional uses of this listing might in- 



nitudes for blended images. Following DargetaL 



( |2010p , who deal with magnitudes for merging sys- 
tems, we use the magnitude rtot of the r-band 
light from both components of each pair, which 



elude correlation with supernova searches, yielding 
complementary reddening information, use of fore- 
ground absorption to distinguish which members 
of galaxy pairs are in the foreground and hence 
constrain their orbital location via redshift differ- 
ences, and use of "wrong- way" redshift differences 
as a probe of peculiar motions. 

Further papers will deal with analysis of the 
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Fig. 4. — Fraction of SDSS galaxies with mea- 
sured redshifts z included in this catalog as a func- 
tion of z. The lower, diverging trace at small z 
shows the result of omitting type B pairs, where 
the smaller background galaxy would be missed at 
larger z. The histogram shows the ^^-distribution 
of all galaxies with spectroscopic redshifts in the 
SDSS (as of DR8), scaled down by a factor lO'^'^ 
to fit on the same axes. Where both redshifts in a 
pair are known, the lower value is used, since the 
closer galaxy is the one under study. The selection 
is most complete at low 2:, exacerbated in this view 
by the incompleteness of SDSS selection of red- 
shifts for very bright galaxies. A plateau appears 
at z = 0.1 - 0.22. Above this, the SDSS parent 
population changes as most galaxies observed were 
color-selected to be luminous red galaxies (LRGs) . 
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Fig. 5. — Component magnitude differences — 
T\ as a function of integrated magnitude Riot of 
sample galaxy pairs, over the well-sampled region 
Tioi = 12.5 — 18. The jagged line is a 51-entry 
running median, showing a gradual decrease for 
fainter pairs in which more nearly equal members 
are preferentially selected. The range r2 — ri < 3.5 
remains well populated over this whole magnitude 
range. 
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Fig. 6. — Distribution of sample galaxy pairs with 
total magnitude Rtot- The Euclidean slope is 
shown for reference in evaluating incompleteness, 
with the caveat that time-dilation effects become 
important in the magnitudes for 2: > 0.1. Subsam- 
ples restricted in component magnitude difference 
Am are shown, with slight offsets for clarity. 
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images of a larger subset of these galaxies from 
optical and ultraviolet images, using these addi- 
tional data to address the spatial distribution and 
reddening law of the dust in the optical and ultra- 
violet ranges. 
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A. Previously known pairs and pairs outside the SDSS 

For convenience of comparison, we list in Table |4] the handful of overlapping pairs noted in other studies 
(including 5 recognized serendipitously in HST images, at redshifts low enough to be included in the catalog). 
Among these, 9 lie within the SDSS sky coverage, all of which were independently found in the Galaxy Zoo 
project. Fig. [7| shows the pairs outside the SDSS DR8 imaging region, collecting color images for ease for 
comparison. These are included in the catalog, and may be recognized by their lack of an SDSS ObjID 
identifier. 
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Fig. 7. — Pairs of overlapping galaxies in the local Universe {z < 0.07) outside the SDSS imaging coverage. 
Where possible, the images are similar in construction to the SDSS color composites. Image sources are as 
given m Table [4j the image of Arp 276 was obtained by Rick Johnson with his 0.25m telescope. These may 
be compared with the SDSS images of both newly found and previously known pairs within the SDSS. 
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Table 1 
Overlapping Galaxy Pairs 



Coordinate name 


SDSS ObjID 


ri 


"^2 


^1 


Z2 




Cross-ID 


SDSS 


J000007.03+081645.1 


587743960499880043 


13.43 


14.84 


0.0387 S 




E 


UGC 12890 


SDSS 


J000046. 97+282407.3 


758874298530726152 


13.42 


14.60 


0.0272 E 


0.0292 W 


F 


UGC 12899 


SDSS 


J000058. 93+285442.1 


758874371533308165 


13.82 


16.17 


0.0230 




Q 


UGC 12901 


SDSS 


J000103. 67+343911.0 


758874373141496140 


14.11 


19.98 


0.0423 SE 




F 




SDSS 


J000140. 21+010531. 2 


587731187814695060 


15.72 


21.66 


0.0611 fg 








SDSS 


J000253. 63+315102.2 


758874300139962790 


18.09 


18.92 






Q 




SDSS 


J000320.46+083707.5 


587743961037078573 


13.78 


15.39 


0.0397 




s 


UGC 10 


SDSS 


J000400. 79+160110.4 


587727223561453822 


15.95 


11092 (S) 




F 






SDSS 


J000401. 52-111027.3 


587727177913008261 


14.75 


11392 (SW) 




S 






SDSS 


J000620. 75-105653. 7 


587727225153978471 


17.73 


18.32 


0.1091 NE 




F 




SDSS 


J000954. 49-050116. 8 


587747122131304642 


16.56 


21.50 






misc 




SDSS 


J001122. 29+062321. 6 


587743795144425534 


12.83 


15.47 


0.0201 W 




S 




SDSS 


J001315. 28+000240. 7 


588015509270822939 


17.03 


17.78 


0.0902 SE 


0.0390 NW 


Q Az 




SDSS 


J001347.12+004612.9 


587731187279134987 


16.56 


21.50 


0.1556 N 


0.1531 S 


E 




SDSS 


J001430. 96+154907.2 


587727180601163792 


16.18 


16.84 


0.0806 S 




S 




SDSS 


J001531. 05-004805. 6 


588015508197343379 


15.15 


17.64 


00685 E 


0.1556 W 


F Az 




SDSS 


J001725. 83-005842. 5 


587731185132109990 


14.09 


20.37 


0.0181 N 




S 




SDSS 


J001747.95-122422.2 


587747073275003048 


14.96 


15.37 






F 




SDSS 


J002305. 73+064506.0 


587744044784156688 


14.25 


17.38 


0.0494 N 




S 




SDSS 


J002315. 77+241310. 9 


587740522933715074 


17.43 


18.5 






X 




SDSS 


J005618. 88+152531. 9 


587724199351025860 


16.01 


16.29 


0.0709 N 


0.0771 S 


s 




SDSS 


J005745. 03-002509. 9 


588015508738801700 


13.51 


14.98 


0.0408 N 


0.0437 S 


E 


UGC 599 


SDSS 


J014145. 47+060151. 9 


587744293354668078 


15.61 


16.90 






Q 




SDSS 


J101800. 73+012116. 3 


587726031700492467 


17.42 


30.0 


0.1265 ctr 




X 




SDSS 


J104524. 96+390949. 9 


587735661016318093 


14.09 


144.3 


0.0254 S 


0.0265 N 


SE 




SDSS 


J104645. 67+114911. 8 


588017703996424196 


10.16 


15.73 


0.0030 SW 




S 


NGC 3368, M96 


SDSS 


J105946. 49+173912.4 


587742865816944739 


14.55 


17.77 


0.0295 W 


0.0294 E 


Q 


Arp 198 


SDSS 


J111605. 47+361410.1 


587739097520799818 


15.24 


16.02 


0.0776 E 


0.0785 W 


SE 




SDSS 


J112116. 82+402043.4 


588017719566467080 


14.31 


16.49 


0.0209 S 


0.1147 S 


F Az 




SDSS 


J112924. 32+415219. 3 


588017721177800860 


15.15 


15.15 






E 




SDSS 


J112943. 31+511415. 7 


587732134846070798 


14.08 


14.29 


0.0335 E 


0.0345 W 


SE 




SDSS 


J113912. 22+553957.8 


587731870706565177 


14.15 


14.42 


0.0623 SW 


0.0615 NE 


SE 




SDSS 


J113947.32+435031.7 


588017625615368313 


16.11 


17.33 


0.1344 W 








SDSS 


J114222. 33+084611. 5 


587732769978777642 


13.91 


13.98 


0.0216 N 


0.0220 S 


F 


IC 720 


SDSS 


J115128. 22+220133.4 


587742061070057546 


13.48 


15.06 


0.0284 W 


0.0257 E 


E 


NGC 3926 


SDSS 


J120802. 07+094557.0 


587732771055337475 


14.52 


15.27 


0.0694 S 




F 


UGC 7114 


SDSS 


J121748. 55+463454. 8 


588298661962973323 


17.84 


20.0 


0.0658 




X 




SDSS 


J123232. 82+635238. 3 


587728676858101872 


14.35 


18.3 


0.0099 N 




B 


UGC 7700 


SDSS 


J123706. 39+234712. 7 


587742188833210435 


15.89 


16.49 


0.0570 S 




Q 


IC 3595 


SDSS 


J124613. 57+171012. 8 


587742774026633226 


14.69 


16.90 


0.0804 NW 




$ 




SDSS 


J125224.50+071053.3 


588017724940878045 


16.68 


17.01 


0.0806 N 




$ 




SDSS 


J130431. 54+303417.7 


587739721901932619 


15.69 


16.98 


0.0622 E 




X 




SDSS 


J131012. 03+373125.4 


587739098604503163 


15.21 


16.82 


0.0429 NE 




B 




SDSS 


J131158. 38+444832.1 


588017627233124523 


14.75 


14.75 






Q 


UGC 8281 


SDSS 


J131206. 27+461146.1 


588298662503972874 


12.9 


17.8 


0.0283 N 




F 


NGC 5021 


SDSS 


J133026. 37+300144.4 


587739709015261237 


14.53 


17.5 


0.0373 NE 




X 


UGC 8497 


SDSS 


J133558. 58+014348.2 


587726031722184751 


14.38 


15.90 


0.0228 S 




Q 




SDSS 


J135316. 39+094017.5 


587736543089000513 


17.33 


18.15 


0.1349 W 




B 




SDSS 


J135507.95+401003.4 


587736585499246631 


14.45 


20.06 


0.0083 NE 




B 


UGC 8841 


SDSS 


J140138. 45-022558. 8 


587729777442488436 


15.08 


17.63 


0.0505 SW 




$ 




SDSS 


J142459. 13-030401.2 


587729776371302561 


14.87 


15.07 


0.516 SE 


0.0513 NW 


S 




SDSS 


J143243. 04+301320. 8 


587739380987592822 


14.52 


16.16 


0.0623 (NE) 








SDSS 


J143243. 90+301329. 2 


587739380987592820 


14.52 


16.16 


0.0623 NE 




Phi 




SDSS 


J144230. 35+222110. 7 


587739809952890969 


14.51 


15.54 


0.0507 NE 




X 


VV 752 


SDSS 


J144750. 76+314553. 3 


587739132421799955 


14.63 


16.58 


0.0557 NE 


0.0454 SW 


SE Az 


IC 4508 


SDSS 


J145013. 12+241919.0 


587739720301871290 


15.26 


15.38 


0.0329 NW 


0.0329 SE 


S 




SDSS 


J161117.22+141531.5 


587739845394366485 


14.01 


14.98 


0.0329 N 




E 




SDSS 


J204719. 06+001914. 8 


587731173842026987 


12.37 


12.63 


0.0140 NW 


0.0127 SE 


F 


NGC 6962 


SDSS 


J220412. 01-083836.0 


587730816822411334 


15.19 


16.86 


0.0640 S 


0.1353 N 


E Az 
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Table 1 — Continued 



Coordinate name 


SDSS ObjID 


n 


r-2 


zi 


Z2 


Type Az 


Cross-ID 
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Table 2: Summary of galaxy-pair subsamples 



Code Description 



Number 



F Face-on spiral and background E/SO 369 

Q Background galaxy edge-on and radial 237 

$ Foreground disk edge-on and background E/SO 156 

X Crossing edge-on disks 200 

SE Spiral/elliptical pairs not otherwise listed 102 

S Two spirals 584 

B Background galaxy has small angular size 181 

E Two E/SO galaxies 59 

R Could also be polar ring 6 

Redshifts indicate not physically associated 218 
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Table 3 

Deeply overlapping galaxy pairs 



Coordinate designation 


SDSS ObjID 




Obj 


Z-2 


Obj 


Type l^z 


Nuc sep" 


Cross-ID 


SDSS 


J011128. 81+261816. 6 


758877156282466365 










S 


4.3 




SDSS 


J020149. 25+001646. 7 


587731513147785324 


0.0434 








R 


< 2 




SDSS 


J072900. 48+430402. 7 


587738066725896649 


0.1885 


abs 


0.1827 


em 


SE 


2.8 






103713.20-274104.0 




0.0095 


fg 


0.0155 


bg 


S 


2.9 


NGC 3314 


SDSS 


J112917.08+353432.6 


587739305286303747 


0.0346 


NW 


0.0345 


SE 


S 


5.8 


NGC 3695 


SDSS 


J113336. 25+560030.9 


587731889502748827 


0.0609 








X 


1.2 




SDSS 


J114512. 95+350510.4 


587739304213938297 


0.0674 








s 


1.0 




SDSS 


J125020. 69+375656.1 


587739098602733640 


0.0351 


ctr 






R 


< 1 


IC 3828 


SDSS 


J131206. 27+461146.0 


588298662503972874 


0.0283 


N 






F 


9.2 


NGC 5021 


SDSS 


J131940. 07+274221. 7 


587741603111632925 


0.0231 








R 


< 1 




SDSS 


J133642. 56+620337.3 


588011219671908454 


0.0311 


fg 


0.01037 


bg 


S 


5.5 




SDSS 


J134134. 82+372625. 7 


588017977826345079 


0.0628 


fg 


0.1713 


bg 


S I^z 


1.2 




SDSS 


J135031. 55+091704.1 


587736542551802124 


0.0653 


N 






Q 


4.7 




SDSS 


J135239. 37+010609.1 


588848901535629486 


0.0716 








R 


< 1 




SDSS 


J150153. 48+353239. 5 


587736586042212586 


0.0496 








R 


< 1.5 




SDSS 


J153322. 82+332933.3 


588017604701257882 


0.0769 


fg 






R 


< 2" 




SDSS 


J155929. 99+094900.6 


587742610274845009 


0.0731 


fg 


0.1025 


bg 


S l^z 


5.2 


CGCG 079-018 


SDSS 


J160737.22+450355.2 


588011101034774582 


0.0441 


W 


0.0442 


E 


S 


5.5 




SDSS 


J162557.25+435743.5 


587729753278578838 


0.0320 


e 






F 


10.1 




SDSS 


J162957.52+403750.5 


587729652348223597 


0.0300 


NW 


0.0289 


SE 


F 


10.2 


NGC 6175 


SDSS 


J163451. 17+481623.7 


587725994646634705 










R 


< 1 




SDSS 


J231444. 84+063821. 3 


587743958884286603 










S 


4.4 
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